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Abstract 

We discuss four-jet production at the LHC. We calculate cross section for both single-parton 
scattering (SPS) using the ALPGEN code as well as for double-parton scattering (DPS) in leading- 
order collinear approach. Our results are compared with experimental data obtained recently by 
the CMS collaboration. We show that the ALPGEN code relatively well describes distributions 
in transverse momenta and rapidity of each of the four jets ordered by their transverse momenta 
(leading, subleading etc.). The SPS mechanism does not explain the distributions at large rapidity 
for the leading jet. The DPS mechanism considerably improves the agreement with the experi¬ 
mental data in this corner of the phase space. In order to enhance the relative DPS contribution 
we propose to impose different cuts. The relative DPS contribution increases when decreasing the 
lower cut on the jet transverse momenta as well as when a low lower cut on the rapidity distance 
between the most remote jets is imposed. We predict very flat distribution in azimuthal angle 
between the most remote jets with low lower cuts on jets transverse momentum. We identify 
phase-space corners where the DPS content is enhanced relatively to the SPS one. 
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I. INTRODUCTION 


The physics of multiparton interactions (MPI), in general, or double parton scattering 
(DPS) for hard processes gained recently a new impulse related to the experiments at 
Large Hadron Collider (LHC) |l-4]. 

The DPS contributions were discused for several reactions to mention double Drell- 
Yan 1^6], creation of bbbb |7], etc. It was predicted |8] and recently found 1^ that the 
double parton scattering leads to very abundant production of (at least) two cc pairs at 
the LHC. Another process where DPS seems to be seen is double //Y production in the 
region of large rapidity separation between the two }Zip's II11I.I12I] . 

We will not discuss in this Letter evident progress in theoretical understanding of the 
DPS or MPI interaction in general, which took place in recent years. In contrast to very 
quantitative and precise studies of many SPS perturbative partonic processes the field 
of MPI and DPS lacks similar precision. So far most of the experimental studies have 
concentrated on the extraction of in so-called factorized ansatz^. In our opinion, at 
present some efforts should be made in order to make the studies of DPS processes really 
quantitative and differential in kinematical variables. In order to make more detailed 
phenomenological studies one needs also more clear cases where the DPS effects could 
be identified without any doubts. 

Four-jet production was traditionally discussed in the context of double parton scatter¬ 
ing. Actually it was a first process where the DPS was claimed to be observed experimen¬ 
tally n4\ . However, in most of the past as well as current analyses the DPS contribution 
to four-jet production is relatively small and single parton scattering (SPS) driven by the 
2 —> 4 partonic processes dominates. To pin down the DPS contribution one has to make 
rather complicated studies comparing Monte Carlo templates and experimental data. In 
such a case the result is not very transparent and one may worry whether the final result 
is not an artifact of an inadequate Monte Carlo generator. A more evident result for DPS 
in four-jet production would clearly provide a new impulse for the MPI community. 

On the theoretical side the DPS effects in four-jet production were discussed in 


Refs. |115^3 i- A first theoretical estimate of SPS four-jet production, including only some 


^ In some cases even experimental extraction of may be subtle 0. 
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partonic subprocesses, and its comparison to DPS contribution was presented in Ref. lIlTI] 
for Tevatron. Some new kinematical variables useful for identification of DPS were pro¬ 


posed in Ref. 11911 . A model dependence on collision energy and minimal transverse en¬ 
ergy oideff was studied in Ref. S. Presence of perturbatiye parton splitting mechanism 


in the context of four-jet production was discussed in Ref. EOII . 

In our recent studies we have shown how big can be the contribution of DPS for (two) 
jets widely separated in rapidity 12111 . Understanding of this contribution is important in 
the context of searching for BFKL effects or in general QCD higher-order effects 12211 . We 
have found that with the present cuts used in the CMS analysis II23I] the DPS contribution 
can be of the order of 10-20%. It could be still somewhat enhanced imposing further cuts 
on transverse momentum of the dijets or azimuthal angle between jets. 

In the present letter we wish to explore exclusive four-jet sample where the situation in 
the context of searching for DPS should be even better. We shall discuss how to maximize 
the DPS contribution by selecting relevant kinematical cuts. 


II. THEORETICAL FORMALISM 


Partonic cross sections used to calculate DPS are only in leading order. The cross sec¬ 
tion for dijet production can be then written as: 




ij—>-kl I 


( 2 . 1 ) 


In our calculations we include all leading-order ij —> kl partonic subprocesses. The 


factor for dijet production is rather small, of the order of 1.1 — 1.3 (see e.g. 1124 l25lh . 


but can be easily incorporated in our calculations. It was shown that already the leading- 


order approach gives resu 


CMS inclusive jet data pllj . 


ts in sufficiently reasonable agreement with recent ATLAS and 


This simplified leading-order approach can be however used easily in calculating DPS 
differential cross sections. The multi-dimensional differential cross section can be written 


as: 


^^DPS(pp 
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where C = 


and partons ], k,l,m = g, u, d, s, u, d, s. The 


2 if hji — hji Akili — ^2^2 

1 if iiji ^ 22/2 V fci/i k2l2 
combinatorial factors include identity of the two subprocesses. Each step of DPS is cal- 


1 (see Eq. (l2.1|l l. 
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2811 arid LHC P29I-I31I1 provide an estimate 


culated in the leading-order approach 
Experimental data from Tevatron 
in the denominator of formula (|2.2|l . In the calculations we have taken in most cases 
= 15 mb. Phenomenological studies of Ceff are summarized e.g. in 13211 with the average 
value (Jgff 15 mb. 

Now we proceed to the SPS production mechanisms of four-jet production. The ele¬ 
mentary cross section for the SPS mechanism has the following generic form: 


d&i 




2 d^PS, 


(2.3) 


where the invariant phase space reads: 

d^ipi d?ip2 d?ip2> 


d^PS = 


d^p4 


{2n)'^5^{pi + p2 + P3 + P4:) ■ (2.4) 


2Ei(27r)3 2E2(27r)3 lE3{lnf 
Above pi, p 2 , P 3 , Pi are four-momenta of outgoing partons (jets). Many possible subpro¬ 
cesses ij —> klmn are possible in general. In some comers of the phase space only some 
processes are really important and others may be safely neglected. 

The hadronic cross section is then given by the integral 


d(T — J dx\dx2 Z fiixi,pl)fj{x2,pp) da, 


ij^klmn 


(2.5) 


ijklmn 


Above fi and /y are parton (gluon, quark, antiquark) distribution function. 

Instead of explicitly using the formulae above we shall use its Monte Carlo version 
as implemented in the event generator ALPGEN [33!]. Weighted events from the gen¬ 
erator will be used to construct distributions presented in the next section. Only light 
quarks/antiquarks are included in our calculation. 


III. FIRST RESULTS 


To start our analysis we wish to check how reliable our SPS tour jet calculation is. In 


Pig. m we confront the results of calcu 


ation with the leading-order code ALPGEN liSSl] 


with recent CMS experimental data 12211 . In this analysis the CMS collaboration imposed 
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different transverse momentum cuts on the leading, subleading, 3'^^ and 4*^ jets (see the 
figure caption). In this calculation we have used an extra K-factor to effectively include 
higher-order effects. How big is the K-factor was discussed recently in 11341] . It was found 
that the next-to-leadtng order contributions are important and lead to i<C < 1. In our cal¬ 
culation here we will use the K-factor found there. We get relatively good description 
of both transverse moentum and pseudorapidity distributions of each of the tour (or¬ 
dered in transverse momentum) jets. Therefore we conclude that the calculation with the 
ALPGEN generator can be a reliable SPS reference point for the DPS effects. 


1 th 




FIG. 1: Transverse momentum (left panel) and rapidity (right panel) distributions of each of 


the four-jets (ordered in transverse momentum) in the four-jet sample together with t 


le CMS 


experimental data 1220. The calculations were performed with the ALPGEN code l33l] . Here 
kinematical cuts relevant for the experiment were applied to allow for a comparison. 


In Fig. 12 we show in addition the calculated contributions of DPS for different values of 
the (Teff parameter. The parameter effectively includes both traditional (independent 
emissions) DPS and effects related to parton splitting l35l. 13611 . In general, the relative 
contribution of both mechanisms may depend on different kinematical variables. Phe¬ 
nomenological exploration of the dependences may help in understanding the interplay 
of the different contributions and shed more light on the underlying mechanism. We ob¬ 
serve that the DPS contribution is rather small (independent oi Ceff parameter) at larger 
transverse momenta and in midrapidities of jets. Their contribution increases when go¬ 
ing to small transverse momenta and large rapidities. We observe a clear improvements 
at large pseudorapidities {\r]\ >3) when the DPS contribution is added to the SPS con- 
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tribution. The presented CMS data were not optimized for search for DPS effects and 
now we wish to explore how to improve the situation, i.e. to enhance sample of the DPS 
events. 




Leading jet t| 


= p |:j ^ 4 jets + X 

(i = 7teV= 


R=0.5 CMS data. 

- SPS; LO(2 
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FIG. 2: Transverse momentum (left column) and rapidity (right column) distributions for the 
leading jet of the four-jet sample. The solid line represents the SPS four-jet contribution, whereas 
the dashed line corresponds to the DPS one. Different values of cr^ff (given in the figure) were 
used in different rows. 
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Having shown that our approach is consistent with existing LHC four-jet data we can 
focus on finding optimal conditions for "observing" the DPS effects. As shown in our 
previous paper on dijets widely separated in rapidity bill the distribution in rapidity 
separation of such jets seems a very good observable for observing the onset of the DPS 
enhancement or even its dominance. In Fig. |3]we show some examples of such distribu¬ 
tions for different cuts on the jet transverse momenta for two collision energies ^/s = 7 
TeV and -^s = 14 TeV obtained with the condition of the four-jet observation. We focus 
only on the distance between the most remote jets; the other two jets are then in between. 
The higher collision energy or the smaller the lower transverse momentum cut the big¬ 
ger the relative DPS contribution. In such cases one can therefore expect a considerable 
deficit when only SPS four jets are included. Such cases would be very useful to "ex¬ 
tract" the (Jeff parameter which for the selected sample does not to be the same as for 
other cases discussed in the literature. Any deviation from the "canonical" value of 15 
mb would therefore shed some new light on the underlying dynamics. For example, a 
two-component model discussed in Refs . 1351.1361] strongly suggests such dependences. 

Another observable discussed in Ref. [12II] is the distribution in transverse momentum 
of the most remote jets. It was argued there that the jets coming from different partonic 
collisions are not correlated and therefore large transverse momenta of such uncorrelated 
jet pairs are possible in general. Here we wish to explore the situation in this respect 
quantitatively for the four-jet sample. In Fig. H] we show an example for one energy and 
one cut on jet transverse momenta. In these calculations we have not included extra cut 
on jet rapidity separation. The DPS and SPS distributions look rather similar for the four- 
jet sample. This shows that an extra cut on the transverse momentum of the pair of the 
most remote jets would almost not help to enhance the DPS contribution. 

Finally in FigJ^we show an example of azimuthal angle distribution between the most 
remote jets. While the distribution for SPS peaks at (p = n, the DPS contribution is very 
flat as the most remote jets come dominantly from different independent and uncorre¬ 
lated (in azimuthal angle) partonic scatterings. Since the DPS dominates we predict very 
flat distribution in relative azimuthal angle. Limiting to small (pjj would further enhance 
the relative amount of DPS without considerable lowering statistics. 

Table I illustrates and summarizes integrated contributions of SPS and DPS for se¬ 
lected transverse momentum cuts imposed on all four jets and on rapidity distance be- 
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FIG. 3: Distribution in rapidity distance of the most remote jets from the four-jet sample for ^/s = 
1 TeV (left column) and \fs = 14 TeV (right column) for different cuts on jet transverse momenta 
(identical for all four jets). 



FIG. 4: Distribution in transverse momentum of the most remote jets for \/s = 14 TeV. 
tween the most remote jets. Already with the lower cut on transverse momentum of 35 
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FIG. 5: Azimuthal correlation between the most remote jets for SPS (dashed line) and DPS (solid 
line) contributions for = 14 TeV. 

GeV (used already in some CMS analyses) the DPS contribution is about 40 % for ^/s = 7 
TeV and about 60 % for ^/s = 14 TeV. Lowering the lower cut on jet transverse momentum 
to 20 GeV gives already 70% and 80% of DPS, respectively. Imposing in addition that the 
distance between the remote jets is bigger than AY > 7 enhances the DPS contribution to 
almost 90%. Therefore imposing such cuts would help to extract fairly precisely the 
parameter from such experimental studies. 


TABLE I: Integrated cross sections in nanobarns for two collision energies and different cuts on 
jet transverse momenta and rapidity distance between the most remote jets. Here, cr^ff — 15 mb 


has been used for calculating the DPS cross section. 


Kinematical cuts: 

\y\ < 4.7 

Vs = 7 TeV 

rrSPS rrDPS DPS 

^ ^ SPS+DPS 

Vs = 14 TeV 

rrSPS rrDPS DPS 

^ ^ 5P5+DP5 

35 <Pt< 100 GeV 

40.55 29.92 42% 

197.74 275.23 58% 

20 < pt < 100 GeV 

1047.37 2 443.77 70% 

4194.11 16 652.39 80% 

20 < Pj < 100 GeV 

AY > 7.0 

18.56 113.26 86% 

151.70 1 194.28 89% 

20 < Pt < 100 GeV 

0 < cpjj < j 

291.68 1221.88 81% 

1 157.15 8 326.19 88% 
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IV. CONCLUSIONS 


In the present paper we have explored how to enhance the relative contribution of 
double-parton scattering for four-jet production. Compared to our previous studies, 
where we focussed on jets with large rapidity separation (important in searches for the 
BFKL dynamics), here we have studied the case of exclusive four-jet production. We have 
shown that already some present data for four-jet production obtained at the LHC by the 
CMS collaboration with relatively small cuts on transverse momenta indicate some evi¬ 
dence of DPS at large pseudorapidities of the leading jet. 

We have shown that imposing a lower cut on transverse momenta and rapidity dis¬ 
tance between the most remote jets improves the situation considerably, i.e.enhances the 
relative contribution of DPS. A dedicated analysis of the DPS effect is possible already 
with the existing data sample at ^/s = 7 TeV. The situation at larger energies, relevant for 
LHC Run 2 should be even better. We predict that azimuthal correlation between jets 
widely separated in rapidity should dissapear in the considered kinematical domain. We 
have found that in some comers of the phase space the DPS contribution can go even 
above 80 %, not necesserily at the expense of lowering the cross section (statistics of ex¬ 
perimental data). 

We have presented the detailed predictions. Once such cross sections are measured, 
one could try to extract the parameter and try to obtain its dependence on kine¬ 
matical variables. Such dependence can be expected due to several reasons (parton- 
parton correlations, hot spots, perturbative parton splitting) and detailed experimental 
data would help us to understand the situation much better. 
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